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I illustrate the phcnomcnological application of Dyson-Schwinger equations to the 
| calculation of meson properties observable at TJNAF. Particular emphasis is given 

Q\ ■ to this framework's ability to unify long-range effects constrained by chiral sym- 

OS ' metry with short-range effects prescribed by perturbation theory, and interpolate 

between them. 
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1. Dressed Quarks 



The Dyson-Schwinger equations (DSEs) 0] provide a nonperturbative approach 
. to studying the continuum formulation of QCD, making accessible phenomena such 

as confinement, dynamical chiral symmetry breaking (DCSB) and bound state 
structure. However, they also provide a generating tool for perturbation theory 
and hence their phenomenological application is tightly constrained at high-energy. 
This is the particular feature of the phenomenological application of DSEs: their 
ability to furnish a unified description of high- and low-energy phenomena in QCD. 
It is elucidated in Refs. ^,|| and here I only illustrate this aspect, using as primary 
^3 ' exemplars the electromagnetic pion form factor and the 7*7r — > 7 transition form 

factor, which are particularly relevant to the TJNAF community. 

A key element in the description of hadronic observables is the dressing of the 
quark propagator, which is described by the quark DSE: 

•i-H . 

><; Sjip)- 1 := i 1 - P A s {p 2 ) + B f { P 2 ) =A f (p 2 )(i 7 -p + M f (p 2 )) (1) 

A 

= Z 2 ( i7 • p + m)™) + Z x J g 2 D, v {p - q )^S f (q)Tl a ( q ,p). (2) 



Here / is a flavour label, D fiV {k) is the dressed-gluon propagator {H), r£ a (g,p) is 
the dressed-quark-gluon vertex, m^ m is the A-dependent current-quark bare mass 

and := J A d 4 g/(27r) 4 represents mnemonically a translationally-invariant regu- 
larisation of the integral, with A the regularisation mass-scale. The quark-gluon- 
vertex and quark wave function renormalisation constants, Z\ and Z2, depend on 
the renormalisation point, £, and the regularisation mass-scale. 

The qualitative features of the solution of Eq. (Q) are known. In QCD the 
chiral limit is defined by rh = 0, where m is the renormalisation-point-independent 
current-quark mass. For rh = there is no mass-like divergence in the perturbative 
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Figure 1: Quark mass function obtained as a solution of Eq. (|^) using D^ v {k) and 
Tl a (q,p) from Ref. |, and m£ >d = 3.7 MeV, m« = 82MeV, m< = 0.58 GeV and 

ml = 3.8 GeV (C = 19 GeV). The indicated solutions of M 2 (p 2 ) = p 2 define the 
Euclidean constituent-quark mass, Mf, which takes the values: M® d = 0.56 GeV, 
Mf = 0.70 GeV, Mf = 1.3 GeV, M b E = 4.6 GeV. 



evaluation of the quark self energy and hence for p 2 > 20 GeV 2 the solution of 
Eq. § is § 
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where 7 m = 12/(33 — 2A^) is the gauge-independent mass anomalous dimension 
and (qq)° is the renormalisation-point-independent vacuum quark condensate. The 
momentum-dependence is a model-independent result. The existence of DCSB 
means that (qq)° ^ 0, however, its actual value depends on the long-range behaviour 
of Dp, u {k) and r° a (q,p), which is modelled in contemporary DSE studies. Requiring 
a good description of light-meson observables necessitates (qq)° « —(0.24 GeV) 3 . 
In contrast, for m/ =/= 0, 



M f (p 2 ) 



2\ large-p 



m f 



(» 



In 



9 2 /A 2 



QCD 



(4) 



An obvious qualitative difference is that, relative to Eq. (Q), the chiral-limit solution 
is l/p 2 -suppressed in the ultraviolet. 

There is some quantitative model-dependence in the p 2 -evolution of Mf(p 2 ) 
into the infrared. However, for any forms of and r£ a that provide an accurate 
description of f n< K an d m^^Ki one obtains profiles like those illustrated in Fig. [TJ. 
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The evolution to coincidence between the chiral-limit and u, d-quark mass functions, 
apparent in this figure, makes clear the transition from the perturbative to the 
nonperturbative domain. The chiral limit mass-function is nonzero only because 
of the nonperturbative DCSB mechanism whereas the u, <i-quark mass function 
is purely perturbative at p 2 > 20GeV 2 , where Eq. (Q) is accurate. The DCSB 
mechanism thus has a significant effect on the propagation characteristics of u, d, s- 
quarks, and this is fundamentally important in QCD with observable consequences. 



2. Bound States 

A meson is a bound state of a dressed-quark and -antiquark, and its internal 
structure is described by a Bethe-Salpeter amplitude obtained as the solution of 

[rp(fc;Q)L = f \xp(q;Q)} sr K™( q ,k;Q), (5) 

where Xp(9>Q) = S{q+)T P (q;Q)S(q^); S{q) = di'Ag(S u {q), S d (q), S s (q), . . .); q+ = 
q + r]Q Q, g_ = q — (1 — T]q) Q, with Q the total momentum of the bound state and 
t]q G [0, 1] the relative- momentum partitioning parameter; and r,. . . ,u represent 
colour-, Dirac- and flavour-matrix indices. For a pseudoscalar meson, such as the 
pion, the solution has the general form 



T P {k-Q) = T F lb 



iEp(k;Q)+j-QF P (k;Q) (6) 



+ j-kk-Q G P (k; Q) + <7 M „ k^Qv H P (k; Q) 

where T is a flavour matrix identifying the meson; e.g., T = i (A 1 + zA 2 ), with 
{V, j = 1 . . . Nj - 1} the Gell-Mann matrices of SU(N f ). 

In Eq. (||), if is the renormalised, fully-amputated, quark-antiquark scattering 
kernel. Important in the successful application of DSEs is that K has a systematic 
skeleton expansion in terms of the elementary dressed-particle Schwinger functions; 
e.g., the dressed-quark and -gluon propagators. The expansion introduced in Ref. ^ 
provides a means of constructing a kernel that, order- by-order in the number of ver- 
tices, ensures the preservation of vector and axial- vector Ward-Takahashi identities; 
i.e., current conservation. Only with such a truncation is an accurate description 
of the light-quark mesons possible. 

In QCD the leptonic decay constant of a pseudoscalar meson is || 

,A 

/ P A1 = trZ 2 / (T p ) T l5l ^ X p(k;Q), (7) 

Jk 

where the trace is over colour, Dirac and flavour indices. Equation (Q) is exact: the 
A-dependence of Z2 ensures that the right-hand-side (r.h.s.) is finite as A — > 00, and 
its £- and gauge-dependence is just that necessary to compensate that of xp(^; Q)- 
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In the chiral limit the axial-vector current is conserved, and employing any 
Ward-Takahashi identity preserving truncation of K one obtains 



f P Ep{k- 0) = B (k 2 ) , F R (k; 0) + 2 f P F P (k; 0) = A Q (k 2 ) , 

G R {k; 0) + 2 f P Gp{k; 0) = 2A' (k 2 ) , H R (k; 0) + 2 /pffp(fc; 0) = , 



where F R , G R and are calculable functions in the dressed axial- vector ver- 
tex, ri* (k;Q). These identities are associated with Goldstone's theorem and in 
fact one can show [|] that when chiral symmetry is dynamically broken: 1) the 
flavour-nonsinglet, pseudoscalar BSE has a massless solution; 2) the Bethe-Salpeter 
amplitude for the massless bound state has a term proportional to 75 alone, with 
the momentum-dependence of Ep(k;0) completely determined by that of Bo(k 2 ), 
in addition to terms proportional to other pseudoscalar Dirac structures that are 
nonzero; and 3) T P (k;Q) is dominated by the pseudoscalar bound state pole for 
Q 2 ~ 0. The converse is also true. Hence, in the chiral limit, the pion is a mass- 
less composite of a quark and an antiquark, each of which has an effective mass 
M E ~ 0.5 GeV. 

For nonzero values of the current-quark mass, whether small or large, instead of 
Eqs. (||) one obtains || 



flm%= -M P {qq)f, (9) 



where M.p = fcff 



M (0 



e.g., for the it: Ai n + — m£ + mj, , and 



{t p ,(t p ) t }" 

~(qq)c = fptvZj (T p ) l5X p(q;Q)- (10) 

Jq 

Equation (^|) is an exact mass formula for flavour non-singlet mesons and I note that 
the r.h.s. does not involve a difference of massive quark propagators: a phenomeno- 
logical assumption often employed. (qq) p in Eq. ([To]) is an "in-hadron" condensate. 
It is gauge-independent and its renormalisation point dependence is exactly that 
required to ensure that the r.h.s. of Eq. (^) is renormalisation point independent. 

For small current-quark masses, rh q ~ 0, Eq. (^J) yields what is commonly called 
the Gell-Mann-Oakes-Renner relation; i.e., m 2 P oc fh q , because 

lim>}f ={qq)l (H) 

However, it also has an important corollary when the current-mass, ttiq, of one or 
both constituents becomes large, predicting || 

m p cx rriQ , (12) 

which follows because (qq) p is TOQ-independent for large-TOQ and fp oc nip 1 ^ 2 . The 
transition from the quadratic to the linear mass-relation occurs at rh q w 2 rh s ^ , 
at which point explicit chiral symmetry breaking overwhelms DCSB. 
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Two other important model-independent results can be obtained [0 from Eq. (||) 
and the systematic construction of K. The scalar functions in Eq. (ra) depend on 
three invariants; e.g., Ep(k; Q) = Ep{k 2 , k-Q, Q 2 )- The zeroth Chebyshev moment 
of these functions; e.g., 



°E P (k 2 ,Q 2 ) :=- [ dx \/l - x 2 E P {k- Q) , k ■ Q := x^fk 2 
n Jq 

are dominant in the description of bound state properties, and 



(13) 



°£p(fc 2 ,Q 2 ) lars « fc " A/ (fc 2 ), (14) 

with °Fp(k 2 , Q 2 ), k 2 °Gp(k 2 ,Q 2 ) and k 2 °Hp(k 2 ,Q 2 ) behaving in precisely the same 
way. Further 

k 20 G P (k 2 ,Q 2 ) larg = " 2 2°F P (k 2 ,Q 2 ) . (15) 

These results determine the asymptotic form of the electromagnetic pion form fac- 
tor. 

3. Electromagnetic Pion Form Factor 
The impulse approximation to the electromagnetic pion form factor is fiof 
{pi+P2)^F v (q 2 ) := A M (pi,p 2 ) (16) 

= 2N c tr D [ r n (k- 1 - P2 )S(k ++ )iri(k ++ ,k + _) 5(fc+_)r ff (fc -;pi) 5(fc__), 

Jk 

where f w (q; — F) T := C _1 r 7r (— g; — P) C with C = 7274, the charge conjugation 
matrix, and k a p := k + ap\/2 + (3q/2, p 2 '■= Pi + No renormalisation constants 
appear explicitly in Eq. ( |16j ) because the renormalised dressed-quark-photon ver- 
tex, TZ, satisfies the vector Ward-Takahashi identity. This also ensures current 
conservation: (p x - p 2 ) l _ l A M (pi,p 2 ) = 0. 

The calculation of F^^q 2 ) is simplified by using a model algebraic parametri- 
sation of the dressed-quark propagator that efficiently characterises the essential 
elements of the solution of the quark-DSE and determines the pion Bethe-Salpeter 
amplitude via Eqs. (||) and fll5|), and an efficacious Ansatz for in which the 
vertex is completely determined by the dressed-quark propagator. The result is de- 
picted in Fig. ^. The current uncertainty in the experimental data at intermediate 
q 2 is apparent in the lower panel, as is the difference between the results calculated 
with or without the pseudovector components: F, G, of the pion Bethe-Salpeter 
amplitude. These components provide the dominant contribution to F- K {q 2 ) at large 
pion energy pc[ because of the multiplicative factors: 7 • Q and 7 • kk • Q, which 
contribute an additional power of q 2 in the numerator of those terms involving F 2 , 
FG and G 2 relative to those proportional to E. Including them one finds 



large-? 2 a{g 2 ) (-(qq)° ^ 2 
I 2 ft 



W) « ^r- : i (17) 
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Figure 2: Upper panel: calculated pion form factor compared with data at small-g 2 . 
Lower panel: the large-q 2 comparison, with the two solid lines showing the range 
of model-dependent uncertainty. In both panels the dashed line fl3|l assumes that 
F '„ = = Gtt = H v , and the data are taken from Refs. £L2l 



i.e., q 2 F 7T (q 2 ) ss const., up to calculable In ^-corrections, in agreement with the 
expectations raised by perturbative QCD. If the pseudovector components of IV 
are neglected, the additional numerator factor of q 2 is missing and one obtains |l3| 
q A Fn(q 2 ) « const.. 

4- 7*7r — > 7 Transition Form Factor 

The impulse approximation to this form factor is 
1 



4tt 2 



(18) 



= tr / S{qi)T 1x {q\-P)S{q 2 )iQ,T1 i {q2,qi2)S{qi2)iQTl{q l 2,qi) , 
J a 
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where Q = (|7 + t 3 )/2 = diag(2/3, — 1/3), and fci, k 2 are the photon momenta 
[on-shell: k\ = = k 2 , 2fci • = -P 2 ], (J is the loop-momentum, and qi := q — ki, 
q 2 := q + k 2 , q := \{qi + q 2 ), <Zi2 := q — k\ + k 2 . The manner in which the "triangle 
anomaly" is recovered with no dependence on model parameters is described in 
Ref. |lOfl , and this same mechanism applies to all anomalous pion and photopion 
processes |l4| . It is a unique feature of the DSE framework. 

Using this expression one can calculate [|l5],[l6) T{k\ , k\-k 2 , fcf ) when one or both 
of the photons is off shell and also determine the asymptotic behaviour analytically. 
The formal character of that derivation is presented in Ref. |l6| but, in neglecting 
essential aspects of renormalisation, it is imprecise. I remedy that here for the 
illustrative case of one photon off-shell: k\ = Q 2 . 

At large-Q 2 , ki ■ k 2 « -Q 2 /2, and in Eq. @ the pion Bethe-Salpeter ampli- 
tude focuses the integration support at q = 0. As a consequence the asymptotic 
behaviour of the integral can be determined using 

S(qi2) ^jr-^H- (h - fca) , ^ » Z lllM , (19) 

which follow from Eq. (j^) and the DSE for the quark-photon vertex, so that 

1 r A 
T^(k 1 ,k 2 )^i—(k 1 ~k 2 ) a trz 2 j ±T 3 xAq;-Ph,*Wv, (20) 

where I have used the Ward identity: Z\ = Z 2 . Hence, the transition form factor 
1 



-^ie^u pa k lp k 2a T(Q Z ) := T^ v (ki, k 2 ) + T^(k 2 ,ki) 
1 ' A 



ie^vpa (h - k 2 ) p § trZ 2 / § r 3 \n(q; ~P)l5la (21) 



— i£p,vpcr k\pk 2a — , (22) 
where the last line follows from Eq. (|7|); i.e., 

T(«V"^"^, (23) 

in agreement with the expectations raised by perturbative QCD. Thus, as with 
Eq. (|l6|), one equation unihes the small- and large-Q 2 results and predicts the 
evolution between them (l3-[l7|. 



5. Epilogue 

I have been necessarily brief. There are many other applications of interest 
to this community, among them the diffractive electroproduction of neutral vector 
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mesons Q , the electromagnetic form factors of their charged states Q and the 
unification of light- and heavy- meson observables || . The most pressing contempo- 
rary challenge relevant to this community is the extension of the framework to the 
calculation of baryon observables, which is underway. 
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